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ABSTRACT 



We use photometric data, from the optical to the 
mm, for a large sample of optically selected radio- 
quiet quasars, at low and high redshifts, to test emis- 
sion models from circum-nuclear dusty torii around 
them. Model parameters, such as dust mass, temper- 
ature distributions and torus sizes are inferred, under 
the assumption that the dust is heated only by the 
central nuclear source. Dependences of best-fit pa- 
rameters on luminosity and redshift are studied and 
the contribution of dust in the host galaxy to the 
observed fluxes is briefly mentioned. Only tentative 
conclusions can be drawn from this yet large dataset 
and detailed modelling, because crucial spectral in- 
formation is missing at wavelengths where the QSO 
dust emission is more prominent, i.e. at A ~ 50 to 600 
fim. Poor information in this wavelength range af- 
fects our understanding not only of the high-redshift 
objects, but also of the low-redshift sample. We then 
expect fundamental contributions by the FIRST mis- 
sion in this field, which is obviously related to the 
more general one of galaxy formation. 
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1. INTRODUCTION 



While disk galaxies appear to form their stellar con- 
tent at moderate rates during most of the Hubble 
time, there are no clear indications yet that the ex- 
pected bright phase of star formation in spheroidal 
galaxies has been detected by optical searches. A 
possible explanation is that galaxies in this phase are 
shrouded in dust, produced in large amounts by the 
early generations of massive stars (e.g. Franceschini 
et al. 1994; De Zotti et al. 1996). 
This idea of an early and quick metal and dust pro- 
duction affecting optical emission at least in bulge- 
dominated massive galaxies may receive some sup- 
port by preliminary evaluations of the dust content in 
high-z quasars and of (very high) metal abundances 
in absorption line systems "associated" to the QSO. 



These have shown that the close environments of 
high-z (both radio-loud and radio-quiet) QSOs were 
the site of quick and efficient star formation and 
metal production (Haehnelt & Rees 1993; Frances- 
chini and Gratton 1997). Indeed, the presence of 
an AGN could even mask the presence of a forming 
spheroid and complicate its investigation, but can 
also be exploited to identify target areas where to 
look for forming structures (galaxies, galaxy clusters 
and groups). 

Sub-mm - mm observations are then starting to shed 
light into these phases. The chance of detection of 
high-z objects at these wavelengths is favoured by 
the strong and positive K-correction implied by the 
steeply rising sub-mm spectra. 

Many high-z sources, most of which associated with 
active nuclei, have been recently detected in the mil- 
limeter (both in the continuum and CO line emission; 
see Andreani et al. 1993; Chini and Kriigel 1994; 
Dunlop et al., 1994; Isaak et al. 1994; McMahon et 
al., 1994; Ivison 1995; Ohta et al, 1996; Omont et 
al., 1996a and 1996b). With the exclusion of radio- 
loud quasars, such radiation has been attributed to 
dust since the observed spectrum, F oc i/ 3-4 , is far 
too steep for any dominant synchrotron self-absorbed 
contribution. 

In this contribution we make use of a large sample of 
optically selected quasars at low and high redshifts, 
observed in the far-IR to mm region, to study dust 
emission in radio-quiet objects and its possible de- 
pendence on redshift. The basic assumption in this 
analysis is that the dust distribution, assumed to 
have simple axial symmetry, either confined to the 
quasar nucleus or covering a sizeable fraction of the 
hosting galaxy, is illuminated by the central nuclear 
source. Obvious evidence in support of this is the 
dominance of the point-like nuclear source in the op- 
tical in both high- and low-redshift objects. 

This effort marks a significant improvement over pre- 
vious studies, as it exploits the whole spectral infor- 
mation from the optical to the mm, the optical pro- 
viding in particular an essential datum on the inten- 
sity of the radiation field illuminating the dust. 
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Figure 1. Observed SEDs of low-z AGNs fitted with the model described in §3 (solid line) with the parameters r " ax , 
6, Ay. (a) 200, 0°, 18 (b)1000, 0°, 40 (c) 1500, 0°, 40 (d) 1500, 0°, 40. Dotted-line curve corresponds to the 
spectrum of a typical spiral galaxy with a dust mass of ~ 10 7 M Q 



2. THE SAMPLE 



3. MODELS OF DUST EMISSION IN QUASARS 



The sample consists of optically selected quasars with 
sub-mm and mm observations. Far-IR data are from 
the IRAS PSC, and for most of the sources co-added 
survey data were provided by IPAC based on the 
SCANPI procedure. Mm observations are partly col- 
lected from the literature and partly taken by the 
authors with the IRAM 30m at Pico Veleta. The 
low-redshift sample consists mainly of PG radio-quiet 
quasars. 



The model spectral energy distributions (SEDs) used 
to fit the observational data were computed with 
the numerical code described by Granato & Danese 
(1994), which solves the radiative transfer equation 
in a circumnuclear dust distribution. This step is re- 
quired since in the torii predicted by unified models 
the dust emission is self-absorbed in the near- and 
mid-IR. The code is quite flexible about the geome- 
try and composition of the dusty medium, the only 
restriction being axial symmetry, and thus allows a 
wide exploration of the parameter space. 

We assume in our exploration that the optical prop- 
erties of dust are the same for all sources, with a 
standard Galactic composition. The inner radius 
ri n of the dust distribution is set by the grain sub- 
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Figure 2. As figure 1 but for high- z AGNs (a) 600, 0°, 40 (b) 800, 0° , 40 (c) 600, 0°, 30 (d) 200, 0° , 18 



limation condition, at T s = 1500 for graphite and 
T s = 1000 for silicates. This translates into the con- 
dition ri n ~ 0.5-^/^46 pc, where L^q is the luminosity 
of the primary optical-UV emission in units of 10 46 
erg s -1 . 

The details of the model, as well as the effect of the 
various free parameters, have been widely discussed 
by Granato & Danese 1994 and in Granato, Danese 
& Franceschini (1997). Here we focus on the simplest 
geometry which is in reasonable agreement with the 
available observations, the "flared disc" , in which the 
scale height of the torus along the z-axis increases 
linearly with the radial distance, and the dust dis- 
tribution within the disc is homogeneous. Model pa- 
rameters are then (i) the outer radius r rnax of the 
dust distribution; (it) the absorption Ay along typi- 
cal obscured directions; and (Hi) the covering factor 
/ = cos(0h) of the torus, where 0/j is the half open- 
ing angle of the dust-free cones. In practice, the free 
parameters in the fit were the former two, the last 
one having been fixed to / = 0.8, a value not incon- 



sistent with requirements of unified models of AGNs 
(Granato et al. 1997). 



Roughly speaking, the first parameter is related to 
the broadness of the IR bump arising from dust re- 
processing, whilst the second controls mainly the 
near-IR slope of the SEDs as observed from obscured 
directions as well as its anisotropy. Therefore Ay is in 
principle constrained by the SEDs of obscured AGNs 
or by testing the anisotropy of mid-IR emission. The 
covering factor is obviously related to the relative 
power reprocessed by the dust in the IR bump and 
the primary optical continuum flux. 
In this simple geometry the total dust mass is 
given by the product between the disc volume 
V = |7rcos0^ and the constant dust density p — 
k jj, T " _ r j , where k v is the absorption coefficient 

per unit dust mass and t v the equatorial optical 
thickness. This gives: 
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Figure 3. Dust masses versus (a) B-band absolute magnitude and (b) redshift. (c) Maximum size of the dust emission 
region as a function of redshift. Filled circles correspond to objects for which the presence of a host galaxy is required to 
improve the fit at A > 60[im. Horizontal lines in (b) show the range of dust masses found in nearby spirals. 
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Unless otherwise stated, all quantities are computed 
assuming Ho=50 km/s/Mpc and a flat universe (qo = 
0.5). 



for comparison. Clearly, while for local objects the 
hosting galaxy can contribute to some extent at the 
longest wavelengths, for high-z quasars the observed 
spectra are far in excess of what would be expected 
from a rcdshiftcd SED of a normal galaxy. This re- 
flects in part the much increased power of the illu- 
minating source, following the well-known evolution 
with z of the primary continuum. 
The overall quality of the fits is remarkably good, if 
we consider that only two parameters (r max and Ay) 
are free to model the spectral shape. 



4. RESULTS 



Figures 1 and 2 show typical spectra of 4 low-z and 4 
high-z sources. Some information on the source and 
model parameters are reported in the figure panels 
and captions. The thick solid curves are the best 
fits SEDs, while the dotted line is the broad-band 
optical-mm spectrum of a typical spiral galaxy, used 



Figure 3a plots the best-fit dust mass, Mj, as a func- 
tion of the B-band absolute magnitude, Mb, show- 
ing the expected trend of increasing Md at increasing 
Mb- A similar trend is observed with the AGN bolo- 
metric luminosity. 

Figure 3b shows the evolution of the AGN versus 
redshift. The filled circles correspond to those ob- 
jects for which the contribution of the host galaxy is 
important. 

Figure 3c plots the best-fit outer radius r out of the 
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dust distribution versus rcdshift. M4 and r out are 
very tightly related in our model, both being precisely 
constrained by the long wavelength (mm or sub-mm) 
flux data. Therefore, the plots in Figs. 3b and 3c 
display quite a similar behaviour. 
In both cases the low-z objects show widely dispersed 
values in the parameters, with typically moderate 
values for the dust mass, Mj <~ 10 5 — 10 7 M© and 
from 10 pc to 1 kpc for r max . For the high-z subsam- 
ple, typical values for Md and r max are larger, with 
Md typically in excess of 10 7 and r max in excess of 1 
kpc. 

Note that the procedure of the dust mass estimation 
is quite a robust one: we believe that uncertainties in 
the model parameters affect its value by less than a 
factor two, and the basic conclusions are not affected. 



5. DISCUSSION 



Some limitations of the present analysis have to be 
considered before drawing any conclusions. 

The first is that the observational SEDs are inter- 
preted assuming that most of the IR-mm flux comes 
from dust illuminated by the AGN itself and that the 
dust in the flared disc is homogeneously distributed 
from the inner sublimation surface up to kpc-scale 
distances. 

For most of the local objects only upper limits to the 
mm fluxes are available. So, it cannot be excluded 
that for some of these we are missing significant con- 
tributions of cold dust from the hosting galaxy at 
A > lOO^m ( most of the mm observations were col- 
lected with large antennas fed with single channel 
bolometers, having therefore a field of view limited 
to a few arcsec). Indeed, Danese et al. (1992) and 
Rowan-Robinson (1995) already suggested that the 
FIR spectrum at A > 60/ito in low-luminosity AGN 
may be contributed by cold dust in the surround- 
ing galaxy situated at large radial distances from the 
AGN and with typical dust masses of 1 to 3 10 7 
M . 

The relative contribution at far-infrared wavelengths 
of the nuclear and extended component will only 
be settled by high resolution mapping with planned 
large mm arrays (e.g. LSA). 

On the other hand, the SEDs of high-z object SEDs 
are less constrained, as only a few spectral data are 
available for each source, due to the fact that IRAS 
was not sensitive enough to detect a significant num- 
ber of them (a couple of exceptions appear in Fig. 
2). In these cases the sub-mm data are the only in- 
formation we have on the circum-nuclear interstellar 
medium, IRAS providing a weak constraint on the 
dust temperature distribution. Comparison of our 
model to the data suggests the existence of a large 
gaseous structure around the nucleus (with 1-5 kpc 
radius), including a dust mass of up to a few 10 8 M , 
and a corresponding gas masses of typically 10 10 to 
10 11 Mq. 

We note, in particular, the large values of the outer 
radius of the dust distribution, which are dictated by 



the combined constraints set by the mm flux, the pri- 
mary optical field intensity and the IRAS upper lim- 
its. Such extended distributions indicate that they 
are not merely circum-nuclear dust torii, as seen in 
local Seyfert galaxies, but have rather scale-lengths 
comparable with those of the bulge of the putative 
host galaxy. 



In principle, we would expect that our estimate of the 
dust mass Md in quasars and the size r max of the dust 
distribution should be unaffected by the power of the 
primary continuum and by its evolution with cosmic 
time. In practice an obvious bias operates in Figure 
3: the more luminous AGNs found at the higher z 
are also likely to be hosted by higher mass galaxies. 
Small mass objects at high-z escape detection and 
the lower right part of the plot cannot be sampled by 
present instruments. 

What is clear is that no high mass objects are found 
in the nearby universe. 



A way to explain the trends observed in Figs. 3b and 
3c would be to assume that we are sampling the same 
class of objects. Evolution in the dust mass content 
in quasars could then be related to the history of 
galaxy evolution, as revealed by optical deep surveys: 
a downsizing process is taking place, such that more 
massive galaxies form at higher redshift, followed by 
a sequence of less and less massive objects forming 
at lower and lower redshift down to the formation of 
dwarfs. 



However, this opens the question of the fate of the 
dust: where has all this dust gone? Was it swept 
away from the quasar? Was it consumed to form 
stars in a late phase of star formation? Note that 
the difference between low and high redshift objects 
would be further enhanced in an open universe: dust 
masses at redshifts larger than 2 for £1 <~ 0.01 are 
higher by a factor between 2 and 5. 



The alternative possibility is that we are observing at 
high- and low-redshifts two distinct populations. En- 
vironmental conditions at high redshifts could favour 
the formation of higher mass black holes, while low-z 
AGNs may be related to the refueling of small dead 
black holes in late type galaxies, where gas is still 
available (Hachnelt & Rees 1993). 



Because of a selection bias (all sources were opti- 
cally selected), no edge-on objects were found in the 
present sample: the nuclear source is not totally 
obscured from the circumnuclear torus. Obscured 
sources, of which IRAS F10214 could represent the 
prototype, can only be found from FIR/sub-mm sur- 
veys such as those expected with the FIRST Satel- 
lite. 

An abundant population of dusty objects could 
heretofore escape detection. Are many early quasars 
obscured by the products of massive galaxy-forming 
starbursts? Such objects might observationally ap- 
pear very much like ultraluminous infrared galaxies, 
so that sub-mm photometric surveys by FIRST rep- 
resent the only way to discover them. 
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